Pseudogap phase in cuprates: oxygen orbital moments instead of circulating currents. 
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Circulating current (CC) loops within the cuprate unit cell are proposed to play a key role in 
the physics of the pseudogap phase. However, main experimental observations motivated by this 
sophisticated proposal and seemingly supporting the CC model can be explained in frames of a 
simple and physically clear microscopic model. We argue that instead of a well-isolated Zhang-Rice 
(ZR) singlet 1 Ai g the ground state of the hole center [CuCi] 5- (cluster analog of Cu 3+ ion) in 
cuprates should be described by a complex 1 Aig- 1 ' 3 B2 g - 1 ' 3 E u multiplet, formed by a competition of 
conventional hybrid Cu 3d-0 2p bi g (a) oc d x 2_ y 2 state and purely oxygen nonbonding O 2p7r states 
with a2 g (jT) and e UXtV {-n) symmetry. In contrast with inactive ZR singlet we arrive at several novel 
competing orbital and spin-orbital order parameters, e.g., Ising-like net orbital magnetic moment, 
orbital toroidal moment, intra-plaquette's staggered order of Ising-like oxygen orbital magnetic 
moments. As a most impressive validation of the non-ZR model we explain fascinating results of 
recent neutron scattering measurements that revealed novel type of magnetic ordering in pseudogap 
phase of several hole-doped cuprates. 



Introduction. The nature of the pseudogap (PG) 
phase to be the most puzzling and anomalous region 
of the phase diagram of the cuprate high-temperature 
superconductors is one of major unsolved problems in 
condensed matter physics. Different theoretical models 
describe the pseudogap state as a precursor of the su- 
perconducting d-wave gap with preformed pairs below 
T* which would acquire phase coherence below T c or a 
phase competing with the superconducting one that ends 
at a quantum critical point, typically inside the supercon- 
ducting dome. 

The order parameter associated with these compet- 
ing phases may involve charge and spin density waves 
or charge currents flowing around the Cu02 square lat- 
tice, such as D-charge density wave or orbital circulating 
currents (see, e.g., Refill for a short overview). In his the- 
ory for cuprates, C. M. Varma— proposes that PG is a 
new state of matter associated with the spontaneous ap- 
pearance of circulating current (CC) loops within Cu02 
unit cell. The current pattern is assumed to disappear 
only at a quantum critical point at a hole doping level of 
x c ~0.19. This Intra-Unit-Cell (IUC) order breaks time 
reversal symmetry (TRS), but preserves lattice transla- 
tion invariance. 

From a theoretical point of view the existence of a 
CC-loop order and the ability of such a q=0 instabil- 
ity to produce a gap in the charge excitation spectrum 
are still highly controversial—^. However, several exper- 
imental observations motivated by this proposal pointed 
to a symmetry breaking in the pseudo-gap phase— ~— and 
provided strong encouragement for models based on CC- 
loop order in copper oxide materials. The TRS violation 
in the PG state of Bi2212 was first inferred from the ob- 
servation of dichroic effect in ARPES measurements^, 
but this measurement has been the subject to a long 
standing controversy (see, e.g., Ref.Hol). 

Seemingly the strongest experimental evidence for 
an orbital-current phase are the observations of an 
unusual translational-symmetry preserving magnetic 



order in underdoped YB^CuaOg+ir, HgBa2Cu04+5, 
La2- 2; Sr 2 ;Cu04, and B^S^CaC^Os+s by spin-polarized 
neutron diffraction 6,8 . However, local probes of mag- 
netism, such as nuclear magnetic resonance (NMR) and 
zero-field muon spin relaxation (ZF-/1SR) have found no 
evidence for the onset of magnetic order at the pseudogap 
temperature T*— ~— . Very recent high-precision ZF-/zSR 
measurements of La2- a; Sr a ;Cu04 in the Sr concentration 
(hole-doping) range 0.13 < x < 0.19— do not support the- 
oretically predicted loop-current phases—, and point to an 
alternative explanation for the unusual magnetic order 
detected by spin-polarized neutron diffraction at lower 
hole doping. Hereafter in this Letter we address such an 
alternative scenario. 

Non-Zhang-Rice model and order parameters. 
The problem of the order parameters in hole doped 
cuprates is closely related to that of ground state of the 
hole centers CuO^ - to be cluster analogues of the Cu 3+ 
ion. The nature of the doped-hole state in the cuprates 
with nominally Cu 2+ ions such as La2Cu04 is a matter of 
great importance in understanding both the mechanism 
leading to the high-temperature superconductivity and 
unconventional normal state behavior of the cuprates. 
In 1988 Zhang and Rice— have proposed that the doped 
hole forms a well isolated local spin and orbital 1 A\ g sin- 
glet state which involves a phase coherent combination 
of the 2pcr orbitals of the four nearest neighbor oxygens 
with the same b\ g symmetry as for a bare Cu 3d x 2_ y 2 
hole. This all assumes that, in the low energy limit, the 
Cu-0 system can be reduced to an effective single orbital, 
or one-band model. 

However, numerous experimental data, in particular, 
recent magnetic neutron scattering findings (see review 
article Ref.Q), suggests the involvement of some other 
physics which introduces low-lying states into the exci- 
tation of the doped-hole state, or competition of con- 
ventional Zhang-Rice (ZR) state with another electron 
removal state. This point was discussed earlier, how- 
ever, mainly as an interplay between ZR singlet 1 A\ g 



and triplet 3 Bi g , formed by additional hole going not into 
big state as in ZR singlet, but into a lff oc d z 2 state—. It 
is worth noting that 3 B\ g state corresponds to a Hund 
3 A 2g term of two-hole e 2 g configuration of an undistorted 
CuOg octahedra. However, later experimental findings 
for very different insulating cuprates and theoretical cal- 
culations have shown that the energy separation between 
the bi g {d x 2_ y 2) and ai g (d z 2) orbitals in CuC>4 plaquettes 
is thought to be of the order of 1.5 eV, i.e. too large for 
quasi-degeneracy and effective vibronic coupling. More 
sophisticated version of the non-ZR states was proposed 
by Varma^, who has proposed that the additional holes 
doped in the Cu02 planes do not hybridize into ZR sin- 
glets, but give rise to circulating currents on O-Cu-0 
triangles. 

On the other hand, cluster model considerations sup- 
ported by numerous experimental data point to a compe- 
tition of conventional hybrid Cu 3d-0 2p b\ g {a) oc d x 2_ y 2 
state with purely oxygen nonbonding O 2p7r states with 
d2g(^) and e UXt y(n) oc p x , y symmetry (see Refs fl6lfl7l and 
references therein) . Accordingly, the ground state of such 
a non-ZR hole CuO|~ center as a cluster analog of Cu 3+ 
ion should be described by a complex 1 Aig- 1 ' 3 B 2g - 13 E u 
valence multiplet with several order parameters such as 
spin and Ising-like orbital magnetic moments, dipole and 
quadrupole electric moments prone to strong vibronic 
coupling, and more subtle hidden order parameters. 

Despite a large body of both theoretical and exper- 
imental argumentation indirectly supporting the exis- 
tence of non-ZR multiplets in cuprates their direct exper- 
imental probing remains to be highly desirable especially 
because there are numerous misleading reports support- 
ing "the stability of simple ZR singlet". For instance, 
the authors of the photoemission studies on CuO and 
F^S^CaC^Os-i— , have reported that they "are able 
to unravel the different spin states in the single-particle 
excitation spectrum of cuprates and show that the top 
of the valence band is of pure singlet character, which 
provides strong support for the existence and stability 
of Zhang-Rice singlets in high-T c cuprates thus justify- 
ing the ansatz of single-band models" . In their opin- 
ion "these states are more stable than the triplet states 
by about 1 eV" . However, in their photoemission stud- 
ies they made use of the Cu 2p 3 / 2 (£3) resonance con- 
dition that allows to detect unambiguously only copper 
photo-hole states, hence they cannot see the purely oxy- 
gen photo-hole a2g and e u states. 

Earlier we have addressed unconventional properties of 
the non-ZR hole center related to the 1 Ai g - 1 ' 3 E u quasi- 
degeneracy (A-E model)-i£. Fig. []] shows the term struc- 
ture of the actual valence A-E multiplet together with 
single-hole basis b\ g {\bi g ) = c d \d x 2_ y 2) + c p \bi g (02p))) 
and e UXtV (\e UXty ) = c v \e UXty (Tr)) + c a \e UXtV (a))) or- 
bitals. The e u orbitals could form two circular current 
p±i-like states, e u ±\ with an Ising-like orbital moment 
(eti±i|Z z|e u ±i) = ±2^^ which is easily prone to be 
quenched by a low-symmetry crystal field with formation 
of two currentless, e.g., p x ^ y -\ike e UXtV states. 
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FIG. 1: (Color online) The term structure of the actual va- 
lence A-E multiplet for a hole CuO^ - center together with 
the single-hole b\ g and e^x.y orbitals. Lower panel illustrates 
the ferromagnetic and toroidal orderings of the oxygen orbital 
magnetic moments within the CuC>4 plaquette. 



Even with neglecting the spin degree of freedom we ar- 
rive at the eight order parameters for the hole [Cu04] 5 ~ 
center including both conventional (two-component in- 
plane electric dipole moment, three-component in-plane 
electric quadrupole moment) and unconventional (Ising- 
like purely oxygen orbital magnetic moment and two- 
component in-plane purely oxygen orbital toroidal mo- 
ment) ones (see Fig.[T]). Thus, the Cu04 plaquette 
with ( A\ g , E u ) valent multiplet forms an unconven- 
tional magneto-electric center characterized by eight in- 
dependent orbital order parameters. Even this simplified 
model predicts broken time-reversal (T) symmetry, two- 
dimensional parity (P), and basic tetragonal (four-fold 
Z4) symmetry. The situation seems to be more involved, 
if we take into account spin degree of freedom, in par- 
ticular, the 1 Ai g - 3 E u singlet-triplet mixing effects. First 
of all, such a center is characterized by a true spin S=l 
moment being gapped, if the ZR singlet A\ g has the 
lowest energy. Strictly speaking, for our two hole con- 
figuration we should introduce two spin operators: net 
spin moment S = §1 + §2 and spin operator V — §1 — §2 
that changes spin multiplicity. It should be noted that 
the y~type order implies an indefinite ground state spin 
multiplicity and at variance with S'-type order is in- 
variably accompanied by an orbital order. The singlet- 
triplet structure of the A-E multiplet implies two novel 
types of the spin-orbital order parameters: spin-dipole 
parameters (Vd x ) and (Vd y ) and spin-toroidal parame- 
ters (VT X ) and (VT y ). Novel ordering does not imply 
independent V-, d- or T-type orders. 

Despite a "fragility" of the orbital e„-currents with 
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regard to a crystal field quenching these can produce 
ferromagnetic-like fluctuations that can explain numer- 
ous manifestations of a weak ferromagnetism in different 
cuprates (see, e.g., Reffl9f) and a remarkable observa- 
tion of a weak magnetic circular dichroism (MCD) in 
YBa2Cu306+a;— . It should be noted that the value of 
MCD effect does not straightforwardly depends on the 
value of orbital, or magnetic orbital moment. Gener- 
ally speaking, the hole doped cuprate could be a system 
with a giant circular magnetooptics if we were able to 
realize the uniform ferromagnetic ordering of the orbital 
e u -currents. It seems likely that the relative concentra- 
tion Xh ~ 10 -4 of circularly polarized e u holes is enough 
to provide the same magnitude of MCD as an applied 
magnetic field of 1 Tesla. It is worth noting that the 
current loop stated, by itself, is incompatible with ferro- 
magnetism and cannot explain the Kerr measurements—. 

It is worth noting that occurrence of both orbital 
toroidal and spin-dipolc order parameters point to 
the hole CuO^ - centers as polar centers with effec- 
tive magneto-electric coupling which can provide ferro- 
electric and magnetoelectric properties for hole-doped 
cuprates 20 . Interestingly, within the 1 A\ gi 1 E u multi- 
plet the electric dipole moment operator can be coupled 
with orbital toroidal and magnetic moments by a remark- 
able magnetoelectric relation 16 i 17 : d x = d me {T y , M z }, 
d y = -d me {T x ,M z }. 

Hereafter we address novel effects related with the 
1 Ai g - 1 ' 3 B 2g quasi-degeneracy (A-B-model). Unconven- 
tional orbital A-B structure of the hole CUO4 hole centers 
with the ground state bl g : 1 Ai g -bi g a 2g (Tr): 1 ' 3 B 2g multi- 
plet (see Fig. [2]) implies several spin, charge, and orbital 
order parameters missed in the simple ZR model. For the 
orbital quasi-doublet 1 A ig - 1 B 2g to be properly described 
one might make use of a pseudo-spin formalism with two 
states 1 Ai g and 1 B 2g attributed to | + ^) and | — 5) states 
of a pseudo-spin s = A, respectively. Then we introduce 
three order parameters: (ct 2 ), (o" x ), and (fry), where &i 
is Pauli matrix. Order parameter (a z ) defines the sym- 
metry conserving charge density fluctuations within the 
CUO4 plaquette: Order parameter (a x ) defines electric 
quadrupole moment of B 2g symmetry localized on four 
oxygen sites: 

Qxy = ^2 ^3/(0 = Qb 29 ■ (1) 

i 

It should be emphasized that the quadrupole moment has 
an electronic orbital origin (see Fig. [5]) and has nothing to 
do with any CUO4 plaquette 's distortions or charge im- 
balance between the density of holes at the oxygen sites. 
It is worth noting that usually a spontaneous imbalance 
between the density of holes at the oxygen sites in the 
unit cell is related to a so-called nematic order. Order 
parameter (a y ) defines an antiferromagnetic (staggered) 
ordering of oxygen orbital moments localized on four oxy- 
gen sites: 

(G z ) = (l lz - l 2z + h z - L) = g L (*y) , (2) 




FIG. 2: (Color online) The term structure of the actual valent 
A-B multiplet for hole CuO^" center together with single-hole 
basis b\ g and a\ g orbitals. Shown are antiferromagnetic (stag- 
gered) ordering of oxygen orbital magnetic moments within 
CuC>4 plaquette (G z -mode) and quadrupole Q xy -moAe. 



where gi w-1.0, if to make use of estimates of the clus- 
ter model—. In other words, maximal value of antiferro- 
magnetic order parameter G z corresponds to a staggered 
order of unexpectedly large oxygen orbital magnetic mo- 
ments m z 0.25 /3 e . In contrast with the net orbital 
moment M z the G z order cannot be easily quenched by 
low-symmetry crystal fields. Fig. [2] shows an illustration 
of the G z order in the CUO4 plaquette. In fact, both 
quadrupole moment Q B 2g and local antiferromagnetic or- 
dering of oxygen orbital moments G z do result from the 
1 Ai g - 1 B2 g mixing effect, in other words, these are a re- 
sult of the symmetry breaking. It should be emphasized 
that the G z order resembles the hotly discussed order of 
circulating currents, proposed by Varma— . 

Two unconventional vectorial order parameters are as- 
sociated with the 1 Ai g - 3 B 2g singlet-triplet mixing effect: 
(VQxy) and (VG Z ). It should be noted that correspond- 
ing orderings do not imply independent (V), (Q X y) or 
(G z ) orders. Moreover, the (VQxy) and (VG Z ) orders 
imply all the mean values (S), (V), (Q xy ), (G z ) for 
CuO^ - center together with their on-site counterparts 

such as (S»), (Qxy(i)), (kz) {i = Cu, Oi, 2 ,3,4) turn into 
zero, at least in first order on the 1 Ai g - 3 B 2g mixing pa- 
rameters. All novel orbital and spin-orbital order param- 
eters appear to be strongly hidden, or the hard-to-detect 
ones and can be revealed only by specific experimental 
technique, for instance, by magnetic polarized neutron 
diffraction. 

Novel orbital modes as seen by neutron diffrac- 
tion. Matrix element of the spin interaction of a neutron 
with a CuO^ - center can be written as follows—: 

47r7? 2 

{SM s T(i\V ppl \S'M sl r'(i') = -r o7 
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(SM s Tii\Y / ^e^\S'M s ,T'fi') ■ (S„ - (e ■ S„) ■ e) , 



(3) 

where sum runs on the two holes (y = 1,2); s;„ is the 
hole spin; e = ^ a unit scattering vector; r electro- 
magnetic electron radius; 7 = —1.913 neutron magnetic 
moment in nuclear Bohr magnetons. Introducing S-V 
representation we make a replacement for the spin mag- 
netic amplitude: 



^ s v e itiT " = is(e tqri + e' iqr2 ) + iv(e 4qri - e lqr2 ) 



where Pi(cos8) is the Legendre polynomial, 9 being the 
angle between q and r. The functions gi(qr) are similar 
to the spherical Bessel functions ji(qr) which appear in 
the expansion of e iq r (see Refl2ll for detail). 

For nonzero orbital matrix elements we obtain after 
some routine algebra 

( 1 A l9 |A(q)| 1 S 29 ) = 



9L , , 
— (cos q x I 



cos q y l)(p x \ [l,/(q-r)} \p y ) 



(10) 



S/s(q)+V/ v (q). 



(4) 



('A^A^E^J = 



In other words, we introduce the S- and y-type spin- 
orbital operators for Cu0 4 ~ centers with corresponding 
order parameters which can be detected by polarized neu- 
trons. As well as spin operator V changes spin multiplic- 
ity, the fv (q) operator changes orbital state. For nonzero 
orbital matrix elements neglecting two-site integrals we 
obtain 

( 1 A lg \f v (q)\ 3 B 2g ) = -^(co S q x l + co S q y l)( Px \e^\ Py ) ; 

(5) 

C-A lg \f v (q\ 3 E ux ^) = T-^^(ctt sin q y , x l(Py,x\e tclr \p x ,y) 



c l(py, x \{lf(<l-r)}\P*,y), (11) 
where (p x \ |i,/(q- r)| \p y ) = -(p y \ |l,/(q- r)| \p x ) and 
(Px\ |/ a ,/(q-r)} \p y ) = 

- 1 {2(g Q (qr)) 2p + (g 2 (qr)) 2p (l - ; (12) 

(Px\ {L,/(q- r)} \Py) = -3i(g 2 (qr)) 2p e x e z ; (13) 



Co- smq x , y l(p Xiy \e lv \p Xiy )) ; 



(Px,y\e lv \Py, x ) = -3(j 2 (qr)) 2p e x e y , 



(6) 



(7) 



{Px,y\e lq - r \p x , y ) = (j (qr)) 2p - -{j 2 {qr)) 2p (t 



x.y 



-■y,x 



), 

where (ji(qr)) 2p is a radial average of spherical Bessel 
function, / = RcuO = a/2, the "-" and "+" signs are 
assigned to matrix elements with E ux and E uyi respec- 
tively. It should be noted that at q = 0: fs(0) = 1, 

Mo) = o. 

Matrix element for the coupling of the neutron spin 
with the electron orbital moment can be written as ([3]), 
if the spin operator Q to replace by an effective q- 
dependent orbital operator as follows— 



A(q) = ^^{U/(q-r w )}e^, 



(9) 



n— 1 v 



where l n „ is the orbital moment operator for i^-hole on 
the n-th oxygen site, n runs over all four oxygen sites in 
the CuC>4~ center, { , } is the anticommutator, 

00 / . \ n 00 

■ r ) = 2 E 7Z o\ = £ il ( 21 + l )aMr)Pi (cose) , 

n=o y ' 1=0 



(Px\ {4,/(q-r)| \p y ) = -M{g 2 (qr)) 



(14) 



The "+" and "-" signs in (fTTj) are assigned to matrix ele- 
ments with E ux and E uy , respectively. Thus, the orbital 
vectorial operator A(q) on the basis of the non-ZR mul- 
tiplet can be replaced by an effective operator as follows: 



A(q) =L G (q)G z + LT(q)T. 



(15) 



In contrast with the spin moment the oxygen orbital 
moment directed perpendicular to the CUO4 plaquette 
in the G z or T XtV modes induces an effective magnetic 
coupling both with z- and x-, and/or y-components of 
the neutron spin. In other words, neutrons see the ef- 
fective orbital magnetic moments to be tilted in c*-q 
plane. Only for q = (go(qr)) 2p = 1, (g 2 {qr)) 2p = 
that is Lg || z, while for nonzero Bragg vectors such 
as q = (OIL) or (10Z) (ga(qr)) 2p and (g 2 {qr)) 2p are of 
a comparable magnitude—, that is \jq can have size- 
able ajy-plane components. For illustration, in Fig[3] we 
present the q-dependence of (go, 2 (qr)) 2p given hydrogenic 
2p functions that clearly supports our message. It is 
worth noting one more that this seeming tilting depends 
both on magnitude and direction of neutron scattering 
vector q. 

We see that the the magnetic polarized neutron diffrac- 
tion measurements provide an unique opportunity to 
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FIG. 3: q-dependence of {go,2{qr))2 P given hydrogenic 2p 
functions (ao is Bohr radius). Arrow point to q- value for 
Bragg vector (101) in YBa2Cu306+i. 

inspect unconventional spin-quadrupole (VQxy), spin- 
dipole (Vd Xt y), locally staggered oxygen orbital (G 2 ), 
and oxygen toroidal orbital (T) orderings, which are de- 
termined by the mixing of the ZR singlet 1 A\ g with non- 
ZR terms 3 B 2g , 3 E U , 1 B 2g , 1 E U , respectively. Exps. (5)- 
(8) and (10)-(14) provide q-dependence of the formfac- 
tors for respective neutron scattering. 

All these novel orbital and spin-orbital orders break 
both the time reversal and tetragonal (Z4) symmetry, 
however, our analysis does show that polarized elas- 
tic neutron scattering measurements performed at the 
Bragg scattering wave vectors can detect two hidden 
modes which preserve translational symmetry of the lat- 
tice, these are the spin-quadrupole (VQ xy )-mode at q = 
(Hi) (see Exps. (5) and (7)) or locally staggered orbital 
G 2 -mode at the Bragg vectors such as q = (OIL) or (10L) 
(see Exps. (10) and (12)-(14)). Namely the latter type 
of the long range magnetic order has been experimen- 
tally observed in the pseudogap phase for three different 
cuprate families, YBCO, Hgl201-£ and Bi2212^. Simi- 
lar short range bidimensional order while occuring at a 
temperature far below T* has been observed in LSCO 
system—. 

Furthermore, to explain the experimental data— £ we 
do not need to engage the spin-orbital coupling, quan- 
tum corrections 2 or orbital currents involving the api- 
cal oxygens— as the measured polarization effects can be 
explained with the locally staggered oxygen orbital mo- 
ments orthogonal to the Cu0 2 planes. It is worth noting 
that the G z -type ordering preserving the translational 
symmetry cannot be detected in the polarized elastic 
neutron scattering measurements performed at the Bragg 
scattering wave vector such as q = (Hi) that does ex- 
plain earlier unsuccessful polarized neutron reports—. 

Oxygen orbital moments must inevitably generate lo- 
cal magnetic fields, first of all it concerns a giant ~ 1 T 
field (given oxygen magnetic moment of ~ 0.1 /is) at the 
oxygen nuclei directed perpendicular to the Cu02 plane. 
However, the 17 NMR data on very different cuprates 
11 do not reveal signatures of static G z type mode. At 
present, there are no published 63,65 Cu or 17 NMR stud- 
ies which give clear results concerning the existence or 



absence of fields of the predicted magnitude in YBCO, 
La-214, Hgl201 or Bi2212. The G z -type orbital magnetic 
order, as any other moment patterns which have reflec- 
tion symmetry across the Cu-O-Cu bonds would gener- 
ate a zero magnetic field on yttrium and barium sites 
in YBa 2 Cu306+(5, YBa 2 Cu 4 08, Y 2 Ba4Cu70i5_(5, thus 
making direct 89 Y and 135 < 137 Ba NQR/NMR methods 
as "silent local probes" despite their pronounced sen- 
sitivity for weak local magnetic fields. This reconciles 
the "non-observance" results obtained by 89 Y NMR in 
superconducting Y 2 Ba 4 Cu70i 5 _5 and 135 ' 137 Ba NQR in 
superconducting YBa 2 Cu40s 12 with neutron scattering 
results—. The ZF-/iSR measurements in YBa 2 Cu306+5 
and La 2 - a; Sr 2 ;Cu04— have also found no evidence for 
the onset of magnetic order at the pseudogap tempera- 
ture T*. The NMR and /iSK experiments clearly rule 
static G z type order out. The failure to detect orbital- 
like magnetic order of the kind observed by spin-polarized 
neutron diffraction surely indicates that the local fields 
are rapidly fluctuating outside the fiSR or NMR time 
window or the order is associated with a small minority 
phase that evolves with hole doping—. 

Summary. Instead of inactive ZR singlet the ground 
state of the hole CuO^ - center should be described by a 
complex 1 Ai g - 1,3 -B 2£ ,- 1,3 E u valent multiplet with several 
unconventional hidden order parameters such as intra- 
plaquette's staggered order of Ising-like oxygen orbital 
magnetic moments or a complex spin-quadrupole order- 
ing. Non-ZR hole centers are believed to be an es- 
sential ingredient of both the hole and electron doped 
cuprates 23 . In particular, these specify many features of 
the pseudogap regime (see, e.g., Refs.0113) which indeed 
manifests clearly seen competing orders. The most con- 
vincing evidence of the non-ZR hole centers in cuprates 
was obtained by experimental findings— which provide 
the first direct evidence of a hidden order incompatible 
with simple ZR singlet. As we argue, firstly one needs to 
reconsider the role of oxygen orbitals, especially the in- 
plane p-orbitals. Indeed, only the competition of the a 
and tt oxygen holes can explain the emergence of oxygen 
orbital moments. It is worth noting that the multiplet 
structure of the ground state makes the hole [CUO4] 5- 
centers in cuprates to be the pseudo-Jahn- Teller cen- 
ters—. Thus, the non-ZR hole centers can be respon- 
sible for different magnetic and lattice effects which are 
often addressed to be signatures of the magnetic (spin- 
fluctuation) or electron-phonon mechanism of the high- 
temperature superconductivity. 

We believe that a large body of puzzling effects gov- 
erned by the non-ZR structure of the hole centers are 
secondary ones and are not of primary importance for 
the high-T c superconductivity. However, we need to un- 
derstand these phenomena to understand high-T c puzzle 
itself. A full exhaustive explanation of all the cuprate 
physics governed by the non-ZR hole centers cannot be 
given at this point, but we propose a consistent picture 
that can be successfully used for the distinctive and com- 
plex description of today and future experimental find- 
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